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Extending the analysis of Q , we obtain a formal expression for the coupling between brane matter 
and the radion in a Randall-Sundrum braneworld. This effective theory is correct to all orders in 
derivatives of the radion in the limit of small brane separation, and, in particular, contains no 
higher than second derivatives. In the case of cosmological symmetry the theory can be obtained 
in closed form and reproduces the five-dimensional behaviour. Perturbations in the tensor and 
scalar sectors are then studied. When the branes are moving, the effective Newtonian constant on 
the brane is shown to depend both on the distance between the branes and on their velocity. In 
the small distance limit, we compute the exact dependence between the four-dimensional and the 
five-dimensional Newtonian constants. 
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I. INTRODUCTION 



Advances in string/M-theory have recently motivated 
the study of new cosmological scenarios for which our 
Universe would be embedded in compactified extra di- 
mensions where one extra dimension could be very large 
relative to the Planck scale. Although the notion of ex- 
tra dimensions is not new [2j, y, La] , braneworld scenarios 
offer a new approach for realistic cosmological models. 
In some of these models, spacetime is effectively five- 
dimensional and gauge and matter fields are confined to 
three-branes while gravity and bulk fields propagate in 
the whole spacetime [S El El- Playing the role of a 
toy model, the Randall Sundrum (RS) scenario is of spe- 
cial interest 9]. In the RS model, the extra dimension 
is compactified on an S1/Z2 orbifold, with two three- 
branes (or boundary branes) at the fixed points of the 
I12 symmetry. In the model no bulk fields are present 
and only gravity propagates in the bulk which is filled 
with a negative cosmological constant. In the low-energy 
limit, an effective four-dimensional theory can be derived 
on the branes 

HHHILH1H1. However, beyond this limit, 
braneworlds models differ remarkably from standard 
four-dimensional models and have some distinguishing 
elements which could either generate cosmological sig- 
natures or provide alternative scenarios to standard cos- 
mology. This has been pointed out in many publications 

E3, EH EE E3, m, EM p3, E3, 113, 13 E3, 12H EE 1^, E3i 

and in particular in [lj where the characteristic fea- 
tures of the model are pointed out in the limit that two 
such three-branes are close to each other. In particu- 
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lar the effective four-dimensional theory was derived in 
the limit when the distance between the branes is much 
smaller than the length scale characteristic for the five- 
dimensional Anti-de Sitter (AdS) bulk. The effective 
four-dimensional Einstein equations are affected by the 
braneworld nature of the model and new terms in the 
Einstein equations contain arbitrary powers of the first 
derivative of the brane distance. In [lj the main results 
have been obtained for the simple case where no matter 
is present on the branes. 

The main purpose of this work is to extend this anal- 
ysis and to derive an effective theory in the presence of 
matter on the branes. At high energy, matter couples to 
gravity in a different way to what is usually expected in a 
standard four-dimensional scenario. In particular gravity 
couples quadratically to the stress-energy tensor of mat- 
ter fields on the brane as well as to the electric part of the 
five-dimensional Weyl ten sor . which encodes information 
about the bulk geometry [23, H3, H3 ■ Consequently 
we expect the covariant theory in presence of matter to be 
genuinely different than normal four-dimensional gravity 
and to bring some interesting insights on the way mat- 
ter might have coupled to gravity at the beginning of a 
hypothetical braneworld Universe. 

Our analysis relies strongly on the assumption that the 
brane separation is small, so that the results would only 
be valid just after or just before a collision. However, it 
is precisely this regime that is of great importance if one 
is to interpret the Big Bang as a brane collision [23, 
I35L l3rj l37| or as a collision of bubbles [3^] • I n particular, 
we may point out 39] where it is shown that bubbles 
collision could lead to a Big Crunch. The authors show 
that, close to the collision, the bubbles could be treated 
as branes. Their collision would lead to a situation where 
the branes are sticking together, creating a spatially-flat 
expanding Universe, where inflation could take place. In 



2 



that model, the collision will be well defined and not lead 
to any five-dimensional singularities. 

In order to study the presence of matter on the branes 
in a model analogue to RS, we first derive, from the five- 
dimensional theory, the exact Friedmann equations on 
the branes for the background. This is done in sectionlTTl 
in the limit where the branes are close together, i.e. either 
just before or just after a brane collision. We then give in 
section IIIII an overview of the effective four-dimensional 
theory in the limit of small brane separation as presented 
in Q and show how the theory can be formally extended 
in order to accommodate the presence of matter on the 
boundary branes. We then check that this theory gives 
a result that agrees perfectly with the five-dimensional 
solution for the background. Having checked the consis- 
tency of this effective theory for the background, we use 
it in order to study the effect of matter perturbations 
about an empty background (i.e. a 'stiff source' approx- 
imation) in section E] For this we consider the branes 
to be empty for the background and introduce matter 
on the positive-tension brane only at the perturbed level. 
We then study with more detail the propagation of tensor 
and scalar perturbations. Although the five-dimensional 
nature of the theory does not affect the way tensor and 
scalar perturbations propagate in a given background, 
the coupling to matter is indeed affected. In particu- 
lar, we show that the effective four-dimensional Newto- 
nian constant depends both on the distance between the 
branes and their rate of separation. We then extend the 
analysis in order to have a better insight of what might 
happen when the small brane separation condition is re- 
laxed. The implications of our results are discussed in 
section fvTI Finally, in appendix iBl we present the tech- 
nical details for the study of scalar perturbations within 
this close-brane effective theory. 



assumed to take their standard fine-tuned values 



A+ = ±- 



k 2 L ' 



(2) 



with any deviations from these absorbed into the matter 
Lagrangians C±. The resulting four-dimensional stress- 
energy tensors on the brane can be written as 



T± = T—9 ± 
^ k 2 L ^ 



r (±) 



In this paper we use the index conventions that Greek 
indices run from to 3 and Latin from 1 to 3, referring 
to the Friedmann-Robertson- Walker (FRW) coordinate 
systems defined below in Q. 

The point of this section is to extract as much infor- 
mation as possible about the dynamics of the system in 
the case of cosmological symmetry in order to obtain a 
result against which the effective theory can be checked. 
Therefore, we both assume the bulk and the brane stress- 
energies t^ v to have the required symmetries. General- 
ising the analysis of |jj , we work again in the stationary 
Birkhoff frame: 



with a 2 (Y) 
n\Y) 



As 2 = AY 2 -n 2 (Y)AT 2 + a 2 (Y)An 2 
C 



-2Y/L 



L 2 a'{Yf 



,-ZY/L 



(3) 



C 



with flat spatial geometry for simplicity. The trajectories 
of the branes are Y — Y±(T) giving the induced line 
elements 



d4 = -(n 2 t -y : 2 )dr 2 +a 2 fc dx 2 



r l)AT 2 + c 2 

j-2 , 2 /. \ i 2 



-At± +a±(i±)dx 2 , 



(4) 



II. FIVE-DIMENSIONAL BACKGROUND 
BEHAVIOUR 



where a±(T) = a (Y = Y ± (T)) and similarly for n±(T). 
The velocities of the branes are completely prescribed by 
the Israel junction conditions |40|: 



We consider a Randall-Sundrum two brane model al- 
lowing the presence of general stress-energies on each 
brane. Specifically, we assume an action of the form 



S 




(1) 



i=± J M 



where the two four-dimensional integrals run over the 
positive- and negative-tension branes A4^ respectively 
and g^ v are the induced metrics. We assume a Z2 sym- 
metry across the branes. 

The five-dimensional bulk is filled with a negative cos- 
mological constant A = —6/k 2 L 2 , where L is the associ- 
ated AdS radius and k 2 the five-dimensional Newtonian 
constant. The tensions X± are, without loss of generality, 



AY ± 2 
F± ( P ±) 



1 



a- 



'^F ± (p d 



k 2 L 
l ±ir p ± 



(±)o 



where p± are the brane energy densities — r 
Hubble parameter on each brane then follows as 



1 



1 



L 2 \F±(p ± ) 
C k 2 
± 3L P± 



L 2 a% 



3G 



Pi- 



(5) 

(6) 
The 

(7) 
(8) 



As in pj we now consider the limit of small brane separa- 
tion by replacing n± and a± with their values no and ao 
at the collision (equivalent to taking the leading order in 
d/L where d is related to the radion, as defined below). 



To this level of approximation the brane position are 
then given by 



Y ± (T) ~ Y oT v ± (T-T ) 



(9) 



v± = n Jl--%F±(p±(T = T )), (10) 
V a o 

where here and subsequently we take ~ to denote the 
leading order in d/L, and we have chosen to consider 
the motion of the branes immediately after a collision 
at T = To and Y — Yq when the branes are moving 
apart. Note firstly that the branes will in general be 
moving with different velocities. Secondly, the limit of 
large energy density p± — > oo corresponds to v± — > no, 
i.e. the limit of null brane velocity. 
The transformation 



T-To 
Y-Yo 



— cosha(y) 
no 

t sinha(y) 
(y — 1) tanh" 1 



(11) 



y tanh 1 



no 



brings the brane loci JSJ to the fixed positions y = 0,1, 
with line element 



ds 2 
d(t) 



d{t) 2 dy 2 -dt 2 + a(y,t) 2 dx 2 
t | tanh" 1 | — | + tanh" 1 [ — 



(12) 



Note as a consistency check that the global coordinate 
t coincides for y = 0, 1 with the proper times t± on the 
two branes (in the small d limit) as defined in e.g. 



dt 
dT 



y=0 



dt+ 
dT 



A generalisation of this metric to 

ds 2 = A{x,y) 2 dy 2 + q^(x, y^^dx", 



(13) 



for branes fixed at y = 0, 1 is the starting point for 
the derivation of the effective theory in the next section. 
There, the proper distance between the two branes is 
measured along a trajectory of constant x M , i.e. it is 
taken to be 




FIG. 1: Comparing coordinate systems: In the near-brane 
limit, the two branes move along the lines Y = Y±(T) (thick 
lines). Lines of constant y (straight) and t (curved) are shown; 
d(t) is the proper distance along a line of constant t between 
the two branes. Note that, if v+ 7^ the values of T at 
either end of this line will be different. The coordinate system 
(y, t) is only defined inside the lines Y — Yq = ±no(T — To) as 
shown 



From (7J and (|10fl it can then be shown that the Hub- 
ble parameter at the time of collision is related to the 
rate of expansion of the fifth dimension with respect to 
proper time t by 



H+(0) 



1 , d(o) 

z tanhAi 



(16) 



— (p+(0)coth d(0) 



p_(0)cosech d(Q) 



Note that this is in general not the same as other defini- 
tions of the radion, more common ones being the distance 
along integral curves of the normal (lines of constant x 
are not in general geodesies in this metric, see Fig. 1) 
or, different again, T_(T) — Y + (T). When the effective 
theory is defined from a moduli space approximation, the 
radion often enters via a ratio of the conformal factors 
on the branes, but this is not meaningful apart from in 
the small- velocity limit. It is however of note that all 
these definitions are proportional in the special case of 
cosmological symmetry and small brane separation. 



d(x) 



dy A{x,yf 



(14) 



In particular, if we choose a specific gauge for which 
A(x, y) is independent of y, the metric (|13|1 is simply 



III. CLOSE BRANES EFFECTIVE THEORY 
DESCRIPTION 

A. Formalism 



ds 2 = d(x) 2 dy 2 + q^{x, y)dx fl dx v . (15) 

As discussed in [l| , it is unclear whether such a gauge may 
be fixed in general, but it can be shown that the resulting 
effective theory is not sensitive to the y dependence of A. 



We work in a frame where the branes are assumed to 
be exactly static at y — 0, 1 with metric (|15|) in order to 
simplify the implementation of the Israel junction con- 
ditions, which would otherwise be difficult. From the 
Gauss equations, the Einstein tensor on a y = const hy- 
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persurface is given by |29j, |4JJ : 



(17) 



The unknown quantity in l|17|) is the electric part of the 
projected Weyl tensor Eg which is traceless, enabling us 
to write the Ricci scalar purely in terms of the extrinsic 
curvature: 



R=-^-K 2 + KiK$. 



(18) 



The Weyl tensor Eg can be expressed in terms of more 
recognisable quantities as 

ES = ~\h K " ~ i D " Dwd ~ K " K + J? 5 *' (19) 

where is the covariant derivative with respect to 
1liv{y), implying from (|17fl that: 



We then have a formal expression for the first deriva- 
tive of the extrinsic curvature in terms of the radion and 
stress-energy: 



dy v 



y=o sinh V O 



cosh VOK {+)t * -K ( -^ 



(25) 



where K £ = Kg(y = 0,1). It is straightforward 
then to obtain the corresponding result for the negative- 
tension brane: 



dy v 



y=i 



d_ 

dy 

d_ 

dy 



Kg 



Kg 



v Kg^ 



E 



„ -'(n-l)! 

+ ( cosh Vd-l] ^-Kg 
y=0 v ' d v 

~ (4> 



v=o 



-sinh VOX 

fd 



sinh V 



coshVO K ( - )tl -K Wtl 



.(26) 



D»D u d+—Kg 
dy 

~Sg + dKKg pi {K 2 KiKf) 



(20) 



where the second line is of higher order in the small dis- 
tance limit d -C L. 

In order to find the derivative of the extrinsic curvature 
4-Kg on the brane, we consider the Taylor expansion: 



oo 1 

Kiv = i) = E - d v n)K " 



(21) 



We expand the n th derivative of the extrinsic curvature 
in powers of d/L, keeping only the leading term: 



gn 



-Kg 



y=o 



= Kg {n) = Kt {n) [ l + 0{- 



and, as shown in Appendix ^ one can obtain the recur- 
rence relation 



where the operator O is defined by 

OZg = [d'"Z a v + d M Z a " - d' a Z" u ] d, 



(22) 



(23) 



giving 



-(2n+l)! " (2n)l 

/C£ (1) -fcoshVo /C£ (0) . (24) 



ri=0 

sinh 



Vo 



B. Einstein equations on the branes 

The next step is to use the Israel junction conditons 
to rewrite the extrinsic curvatures of the two branes in 
terms of the stress-energy tensors and the tensions: 







V 


L 


K t-)n = 






L 




T (±)M . 



2 
2 

_ I t (±)jp 
3 " 



(27) 
(28) 
(29) 



This gives us both the value of the Ricci scalar on the 
branes JTSJl 



R{ ±) = ^ T i± ) + ^U±)2_^) r ^)P 



(30) 



and, substituting (t25|) (or ffify) into pU)l. the effective 
Einstein equations: 



G 



}-D^D {± \d± -L|9d|tanhl^^ (31) 

1 1 A n Ifldl |0d|\ au ,_ , 
— t^-t tanh - — - + tan - — - d^dd v d 
Ld\dd\ V 2 2 / 



± 



.(±> 



+ 2^ 



2k 2 



(±)m 



-(±) 



-Sg, 



K 4 n (±)M 



where 



O cothVOf 



cosechv Of 



(32) 
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and 



with the equation of motion for d: 



n 



(±)M _ 1 (±)a (±)/5™ _ 1 (±)_(±)^ _ J_ T (±)2^ 

(33) 

From the tracelessness of Effi we obtain two equivalent 
equations of motion for the radion, 



□Wd ~ T Mf 3 tanh^-tan^)-^34) 

1j \ Z Z / Z 



i 



_ f (±W(±V 

4 "A 4 



The terms in square brackets in l|31|) and (|34() will turn 
out to be of higher order as d — > and so should not 
strictly be included. However, for exact cosmological 
symmetry, they are the only higher order terms and we 
shall keep them for the time being. Later on they shall 
be discarded. 

The conservation of the stress-energy tensor on both 
branes follows directly from the Codacci equation [2^, 
El: 



D^S - D V K 



0, 



which, evaluated on the branes implies 



D (+) r (+)M 



£><-V (-) " 



0. 



(35) 



(36) 



C. Low-energy limit 

As a first consistency check of this close-brane theory, 
we consider its low-energy limit and compare it with the 
effective four-dimensional low-energy theory [Tol fill Il2l 
Il3l | for small brane separation. In that common limit, 
both theories should agree. 

In the low-energy limit, the magnitude of the stress- 
energy tensor on the brane is small compared to the brane 
tension. Any quadratic term k, 4 t^ 2 is negligible com- 
pared to j-t^, so that IL^' 1 * may be dropped in l|31|) 
and, from <|3T)jl . the Ricci tensor on the brane is: 



R (±) = " C 



(37) 



Furthermore in the low-energy limit, the branes are mov- 
ing slowly, \dd\ <C 1, to linear order in |<9d|, we have: 



A 



(±)m 



f (+)/* 



.*(-)/*_ 



(38) 



The effective Einstein equation on the brane at low en- 
ergy is therefore 

g (±) : ~ i#)^(±)d+^(fW:+f(-':) (39) 



□Wd - ±~ (ddf + ^ (V+) + 



.(-) 



(40) 



We can therefore write (|39|l in the more common form: 



G 



(±)m 



-f 

2d V 



(41) 



which is precisely the result we get from the effective low- 
energy theory in the close brane limit 0, ITTl Il2l IT3I |20| . 



IV. COSMOLOGICAL SYMMETRY 

In the most general case, the coupling of the radion to 
matter on the branes given by (|32[) is intractable. How- 
ever, we are concerned here with the case of cosmological 
symmetry as a check on the validity of the theory. We 
take as our metric and notice that 



O 



A 

BS] 



-A 
B5\ 



(42) 



We can then obtain the coupling tensors in closed 
form: 



A^l = —dcotd ( +p± 



-d cosec d ( —p T + p T 



(43) 



A {±) \ = I (p±coth d + ^cosech ctj 5). (44) 

The resulting equations of motion follow from l(3*T |l .l|3*l ^ 
and (jnni : 



-dH ± ±-t a nh- + —A {±} 
L 2 6 



(45) 



2 t\i 1 A O 

± 3L P± + 18* P± 



d + 3H + d 



3G 



did d 

— 3 tanh tan — 

L \ 2 2 



K 2 d /p+ 



(f--)4( 



Pi 



L V 3 
-3H± (p± +p± 



P± (P± + 3 P± )(46) 
(47) 



H ' A 
2 \ ! 



2 1 2 



3 P+P " 



(48) 



6 



Equations l|48|l and (|46|l together imply 
C 



Hi 



(49) 



where C is an integration constant which, at this order, 
can be identified with the bulk parameter C via Q ■ The 
system is now manifestly finite as d — * 0. Note that, 
apart from the presence of quadratic terms, l|46|l is the 
same result as that obtained from the moduli space ap- 
proximation and is, in fact, exact (d decouples as a con- 
sequence of the simplicity of the Weyl tensor for an AdS- 
Schwarzschild bulk). However, the additional informa- 
tion from l|47l) gives additional information not obtain- 
able from a low-energy effective theory. Since H + takes 
a finite value at the collision, the coefficient of d^ 1 in 145|) 
must vanish at d = 0; this implies then that 



H + (0) = ^tanhM 



(50) 



K 2 / 

+ ~6~ \ p+ ^ coth + /°-(°) cosecn d (°) 
in agreement with the exact result 



V. EFFECTIVE THEORY FOR 
PERTURBATIONS 

More interesting is the study of cosmological pertur- 
bations, for which a relatively straightforward solution 
of the above system is also available. We shall give a 
few examples here and point out some interesting fea- 
tures. Throughout we work only with the positive- 
tension brane, assuming the negative-tension brane to 
be empty, and drop the ± signs. Also, we shall assume 
for simplicity that matter on the positive-tension brane 
is only introduced at the perturbed level, i.e. the back- 
ground solution a(t), d(t) is that obtained from (|46|l and 
(|47|l in the absence of matter. We therefore have 



7#(x,t) = 



d(x, t) 



d(t) 
9^ 



Srjf 

6d\x, t) 



where is the usual flat FRW metric with scale factor 
a(t). However, in the following we will set Sd = 0, either 
because we are considering tensor perturbations only or 
because we choose to work in such a gauge. Hence we 
shall assume that d takes its background value. 



A. Tensor Perturbations 

As the simplest starting point we consider perturba- 
tions using the above formalism and we choose to work 
in conformal time. We take the metric to be 



ds 2 



a(r)) 2 (-di] 2 + (Sij + hij) dx i dx j ) 



(51) 



with the usual transverse traceless conditions 



3,1 







on the perturbation, spatial indices being raised by 5 1J . 
The resulting Ricci tensor perturbation is then 



SR 



SR 0l = 0, 



SRij = (W + 2H 2 ) h tJ + WJ t . - -Dh i:j , (52) 

where TL = a' /a = aH, primes denote differentiation 
with respect to conformal time rj and □ = —d 2 + did 1 
is the Minkowski space wave operator. We assume that 
these gravity waves are sourced by tensor matter at the 
perturbative level, i.e. 



Srt 





r| 



0. 



The perturbed Klein Gordon equation for tensor matter 
just reduces to S (Dd) ~ 0, so it is consistent to set the 
scalar perturbation Sd to zero, i.e. to study purely tensor 
fluctuations. In this case, the equation of motion for the 
perturbations follows from (f5T|) : 



3<5 (DiDj d) + — tanh ~ a 2 h i:j (53) 

IjCL 



d 
2 



H — -V Ocoth VO Stu, 
2d 

where we have now dropped the sub-dominant matter 
terms. It is straightforward to obtain 



i) (D,D , ill = - [ l -h' l3 



G?COthd Tj. 



Hhj I d' 



O coth V O Sf, 



Equations (|45|l and l|46(l then imply the relation for the 
background Hubble parameter 



H' + H 2 = 0, H 2 



d! o d d 
-— Ti + a 2 — tanh— . 
d L L 



Putting this all together we obtain, to leading order in d: 



□ ha = 



d dr] 



hit == — k — cothd 8t4. 



(54) 



The same calculation repeated subject to the low- 
energy approximation, not assuming small d, is straight- 
forward. Since the matter is traceless, the standard equa- 
tions at low energy [HEl [H Eg give 

Rij « ^ (2LD[ +) d ] d-g lJ dd 2 + 2d i dd j d)+^r^ ) , 

where -0 = 1 — exp (— 2g?/L). Perturbing this gives 
(1-0) d r 



n 







2k 2 
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using the equations 



Hi + H 2 = 0, H 2 



, (1-0) 
Lip 



Hd! 



d' 2 (1 - 1/>) 
L 2 V : 



for the background. The small-<i limit of (|55Jl . where 
ip ~ 2d/L, is then 



n; d \' 



The operator □ defined in l|54|l is therefore the same 
as one would find in the low-energy theory. The dif- 
ference lies in the source term; in the high-energy theory, 
the effective four-dimensional Newton constant on the 
positive-tension brane is related to the five-dimensional 
one by 



H 



(+)2 
id 



■ cothti k 2 , 



whereas the low-energy result has 

„2 



„( + )2 
l 4d 



K 

d 



(56) 



(57) 



As is the case in the low-energy effective theory, the cou- 
pling to matter is different for the background as it is 
for the perturbations - for the background, the coupling 
can be identified from (JSJ or (|49|l as k 2 /L, as opposed 
to (|56l) . When cither the branes are stabilised, and d is 
not treated as a dynamical variable, d ~ do = const, or 
the velocity is small d< 1 (which is the case in the low- 
energy limit), it is easy to see that l (5o] l and (jETjl agree. 
However, for arbitrary brane velocities, when the radion 
is not stabilised, the exact result for small d is given by 
(|56[) . As expected, the effective Newton constant picks 
up a dependence on d, as it does in the low-energy the- 
ory, but more unexpected, it also contains some degree 
of freedom: the brane separation velocity. Whilst this is 
not expected to be relevant today, since one would as- 
sume the radion is stabilised in the present Universe, it 
would be extremely important near the brane collision. 
As discussed in section^ d would be approximately con- 
stant, d ~ v, leading to 



„(+)2 



cothw 
t 



(58) 



where the coefficient cothv could take any value greater 
than 1 depending on the matter content of the branes. 



B. Scalar Perturbations 

We now consider scalar metric perturbations on the 
brane sourced by a perfect fluid at the perturbative level 
(again, the background geometry is taken to be empty). 
We choose to work in a gauge where 5d — 0, i.e. to 



evaluate the perturbations on hypersurfaces of constant 
d, in which the metric perturbation can be taken as 

ds 2 = a(r]) 2 ((-l + 2<f>)dr 1 2 +4E^dr ] dx i (59) 



The calculations are not nearly so straightforward as for 
tensors and have therefore been relegated to Appendix 
IbI The result is the following equation of motion for the 
curvature perturbation 



□ 



d/_d_ 

d drj 



* = -a 1 



6 d 



coth dSp, (60) 



giving rise to the same relation between the four- 
dimensional Newtonian constant and the five- 
dimensional as in H56|) . Here again we may check 
that, apart from the modification of the effective Newto- 
nian constant on the brane, the perturbations propagate 
in the given background exactly the same way as they 
would if the theory were genuinely four-dimensional. 
This is a very important result for the propagation of 
scalar perturbations if they are to generate the observed 
large-scale structure. The five-dimensional nature of the 
theory does affect the background behaviour but on this 
background the perturbations behave exactly the same 
way as they would in the four-dimensional theory. 

This result is of course only true in the close-brane 
limit, for which the theory contains no higher than second 
derivatives, only powers of first derivatives. When the 
branes are no longer very close to each other, the theory 
will become higher-dimensional (in particular the theory 
becomes non-local in the one brane limit). The presence 
of these higher-derivative corrections (not expressible as 
powers of first derivatives) is expected to modify the way 
perturbations propagate in a given background, mainly 
because extra Cauchy data would need to be specified, 
making the perturbations non adiabatic pol |. However if 
we consider a scenario for which the large-scale structure 
is generated just after the brane collision, the mechanism 
for the production of the scalar perturbations will be very 
similar to the standard four-dimensional one. 



Relation between the four- and five-dimensional 
Newtonian constant 



The relation l|56|l between 



„(+)2 
v 4d 



and the five- 
dimensional constant k 2 is formally only valid for small 
distance between the branes. However if we consider the 
analysis of , we may have some insights of what will 
happen if we had not stopped the expansion to leading 
order in d. Here, terms of the form dd and more gener- 
ally any term of the form d n d^ n+1 ^ have been considered 
as negligible in comparison to d and therefore only the 
terms of the form d n have been considered in the ex- 
pansion. In a more general case, when the branes are 
not assumed to be very close to each other, any term 
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of the form d n S n+1 ^ should be considered and would 
affect the relation between the four-dimensional Newto- 
nian constant and the five-dimensional one. For moving 
branes, we therefore expect the relation between n 4d 
and k 2 to be: 



,(+)2 



„(+)2 

Md 



L 



n 



L 



(61) 



The relation is therefore a functional of d: f2 [d(t)] has an 
infinite number of independent degree of freedom. 

In the low-energy limit, or when the radion is sta- 
bilised, d ~ do = const, the exact expression of fi is 
421: 



n -> n [d(t) = d ] 



oda/L 



2 sinh do/L 



(62) 



For close branes, another limit is now known: when d <C 
L, 



n [d < l] 



■ cothd. 



(63) 



But in a general case, 51 (and therefore «V/ ) is ex- 
pected to be a completely dynamical degree of freedom. 
For the present Universe the radion is supposed to be 
stabilised, but in early-Universe cosmology, the effective 
four-dimensional Newton constant could be very differ- 
ent from its present value. It might therefore be inter- 
esting to understand what the constraints on such time- 
variation of the Newtonian constant would be and how 
it would constrain the brane velocity 0, 13 > or whether 
such a time variation could act as a signature for the 
presence of extra dimensions. 



VI. CONCLUSIONS 

In the first part of this work, we derived the exact be- 
haviour of FRW branes in the presence of matter. The 
characteristic features come from the presence of the p 2 
terms in the Friedmann equation and from the 'dark en- 
ergy' Weyl term. In the limit of close separation we re- 
lated the contribution of the Weyl term to the expansion 
of the fifth dimension. We then used this result to test 
the close-brane effective theory that was first derived in 
but now with matter introduced on the branes. For 
this we have shown how matter can be included using 
a formal sum of operators acting on the stress-energy 
tensors for matter fields on both branes. In the gen- 
eral case the action of this sum of operators on an ar- 
bitrary stress-energy tensor would not be available in 
closed form, although one could in principle proceed per- 
turbatively. When a specific scenario is chosen, however, 
one can make considerable analytical progress. Assum- 
ing cosmological symmetry, the action of the operators 
on the stress-energy tensor is remarkably simple and the 
sum can be evaluated analytically. We then compared 



the result with the exact five-dimensional result in the 
limit of small brane separation. As expected both results 
agree perfectly. Furthermore we have checked that, in the 
low-energy limit, our close-brane effective theory agrees 
perfectly with the effective four-dimensional low-energy 
theory, giving another consistency check. 

We then used this close-brane effective theory in or- 
der to understand the way matter couples to gravity at 
the perturbed level. In order to do so, we considered a 
scenario in the stiff source approximation for which the 
background is supposed to be unaffected by the presence 
of matter and considered the production of curvature and 
tensor perturbation sourced by the presence of matter 
fields on the brane. Although the five-dimensional na- 
ture of the theory does affect the background behaviour, 
we have shown that for a given background the pertur- 
bations propagate the same way as they would in a stan- 
dard four-dimensional theory This is only true in the 
limit of small brane separation and is not expected to be 
valid outside this regime. However, since the large-scale 
structure of the Universe might have been produced in 
a period for which the branes could have been close to- 
gether (for instance just after a brane collision initiating 
the Big Bang), this regime is of special interest. The 
fact that the perturbations behave the same way, for a 
given background, as they would in a four-dimensional 
theory is a remarkable result for the production of the 
large scale structure which could be almost unaffected 
by the presence of the fifth dimension. On the other 
hand, the relation between the four-dimensional Newto- 
nian constant and the five-dimensional one is however 
affected by the expansion of the fifth-dimension. It has 
been shown in the literature 9, 42] that four-dimensional 
Newtonian constant was dependant on the distance be- 
tween the branes, giving a possible explanation of the 
hierarchy problem. In this paper we show that the four- 
dimensional Newtonian constant also has some depen- 
dence on the brane velocity which we computed exactly 
in the small-distance limit, which might be able to pro- 
vide an observational signature for the presence of extra 
dimensions. Outside the small d regime, we expect the 
four-dimensional Newtonian constant to depend on the 
five-dimensional one not only through the brane sepa- 
ration velocity d but also on higher derivatives of the 
distance between the branes d^ n \ making the require- 
ment for moduli stabilisation even more fundamental for 
any realistic cosmological setup within braneworld cos- 
mology. 
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APPENDIX A: LEADING ORDER DERIVATIVE 
OF THE EXTRINSIC CURVATURE 

In this appendix, we shall derive an expression for the 
derivative of the extrinsic curvature on he branes. We 
will not be able to calculate these quantities exactly, but 
will be able to obtain a relatively simple expression for its 
leading-order contribution. We will focus on the positive- 
tension brane first, and our starting point shall be the 
Taylor series 



OC 



(Al) 



where we are defining 



d'' 



" ~ dy n v 



y=o 



We are interested only in the leading order contribution 

toi^ (n) , 

X7f(») =K.^ n] (l + 0(d/L)), 
and the aim of this section will be to establish that 

K,^ n) = OK.^ n - 2 \ (A2) 
where the operator O is defined by 

OZ£ = [d'^ZI + dvZ^-d^Z^dc. 

This implies that JCu is of the same order as /C^" ^ , 
and will allow us to produce a simple, albeit formal, ex- 
pression for this sum, which will be the starting point 
for writing down the small-d effective theory in the next 
section. We will proceed by induction, and throughout 
make the following assumptions about the order of terms: 

• dfj,d v d ~ d v d ~ d° 

• t !*, Da V are at worst as divergent as the 
geometry 

. E i+)li ~d- 1 . 

The assumption on the order of magnitude of the mat- 
ter terms is reasonable, since the matter introduced on 
the brane is expected to scale as the scale factor for the 
background and as the curvature perturbation for gen- 
eral perturbations. Since the curvature perturbations is 
in general expected to diverge logarithmically at the col- 
lision, we can hence assume that t' + ^, is, at worst, loga- 
rithmically divergent in d. This implies that the extrinsic 
curvature on the brane is itself at worse logarithmically 
divergent in d. Similarly, we know that E (+)t * ~ d" 1 for 
the low-energy theory. Although we have argued that 
at high energy the moduli space approximation does not 
give the exact expression for the Weyl tensor, we have 



seen that (at least for the background) the behaviour is 
the same, differing only in corrections at higher order in 
the velocity. In particular Eff should go as d -1 at high 
energies as well (we will see later that this is indeed the 
case). From (I19|) we have 



= -dE* (y = 0) - D»d v d\ y=Q 

L v 4 a 



(A3) 



which, from the above assumptions, gives us 

/C£ (1) ~ dP. (A4) 

For the second derivative of the extrinsic curvature, i.e. 
for ii = 2, we need expressions for the derivatives of the 
Weyl tensor and the Christoffel symbols. It is straight- 
forward to show that 

T% = D^{dK°) + D v {dK«) - D a (dK,„), (A5) 

and the derivative of the Weyl tensor is 

£f =d{2K^E^-^KEi;-^E^ + C^K^ (A6) 

2K«K a pk° - Ik^k^k^ - x - qiiV k a!i k< 3 p k a p 

-Ln a [d 2 D»K av + d 2 D u K£ - 2d 2 D a K»] , 

where kg — K£ — \KS£. On the brane, from the Israel 
matching conditions, the trace of the extrinsic curvature 
is K ^ t, hence k% <~ r also. So the cubic terms in k 
will be of higher order than the KE terms, as will the 
CK term. The leading terms will, in fact, just be the 
first three, giving 



(A7) 



On the brane, 



D 



= {d a d)K (+) » + dD^K {+ ^ 



(d a d)K^, 



(A8) 



the second term being subdominant from the assumption 
that D { a +) K {+) ^ cx L»i +) f (+) ^ is of higher order than cT 1 . 
The derivative of the Christoffel symbol will similarly be 
of the same order as the extrinsic curvature on the brane: 

r°'(0) ~ (d»K {+) : + d „K {+) ; - d' a K$) . (A9) 
Taking the derivative of H±9|l gives 

Kf{y) = -dE£' +2dK^D p d v d+q^T^d a d 

-d d y (K£K«) , (A10) 
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in the bulk. Evaluated on the brane using i|A7(l and l|A9[l . 
the dominant term (of order d°) is the one containing the 
derivative of the Christoffel symbol 



(All) 



Since we have shown in (|X9|) that P^(0) ~ dj 3 K {+) ", on 
the brane, the second derivative of the extrinsic curvature 
is hence of the same order as the extrinsic curvature itself 
K^ 

Using (|A9|> . we have proved the result for n = 2: 



OK 



d^K i+] * + d u K {+)a " 
(+)m 



(A12) 



The second derivative of the Christoffel symbol follows 
from HASH: 



r M V'(y) = D^dK^ + D^dK^-D^dK 1 ^) 

+d (r%KZ + T%K> - 2T^K«) (A13) 



1 iSv-^p 



where, recall, = dC n K^ is a tensor, hence the use 
of covariant derivatives. These D[dK') terms are of or- 
der d°, whilst the others are all of higher order, when 
evaluated on the brane. The leading term is 

r^"(0) ~ d^dK a u '(Q) + d„dK«'(Q)-d a dK^(Q) 

~ d°. (A14) 

Substituting (|A6|I into (|A10|) gives a complicated second- 
order differential equation for K£. Taking repeated y- 
derivatives of this equation would be impractical, but to 
start with all we want to do is to work to leading order. 
We will first identify which term is dominant, before ac- 
tually evaluating it. We therefore drop all indices and 
numerical factors for the time being, writing q for the 
metric (with indices in any position), K for Kj} and d 
for <9 M . For example, q' oc dK^ v = dq^K® and so we 
would write q' = dqK. The equation for K" can then be 
written symbolically as 

K" = d (d 2 + Td + dT + T 2 ) q' + dK 3 + dK (A15) 
+dKK' + qdd(T' + dTK + ddK) + d 2 dq' 1 

and, from (|Allfl . we already know that the dominant 
term is qddV ~ K. We know that K im \0) and T im \0) 
are all of order d° or K(Q) for m = 0,1,2. Recalling 
that the extrinsic curvature on the brane is at worse log- 
arithmically divergent in d, terms of the form dK(0) will 
hence be negligible compared to terms of order d° (and 
of course compared to terms of order K(0)). Compared 
to the dr 1 divergence, K(0) is hence still negligible. In 
what follows, terms of order d° (such as ^'(0), T(0) and 



r"(0)) and terms of order K(0) (such as K(0), K"(0) and 
r'(0)) can hence be treated in a similar way. Since they 
are all at worse going as K (0) , we shall use in what follows 
the notation K^(0) ~ r( m )(0) - K(Q) for m = 0,1,2. 
We shall hence take as the inductive hypothesis that this 
result is true for all m < n. In particular, 



T (m) (0) - K(0) V0<m<n 
(0) - dK(0) V < m < n 



Jm+l) 



g(0) ~ d°. 
Writing I = n — 1 , we have 

K (n+1) = Q(l) K „ 



(A16) 



d (d 2 + Td + dT + T 2 ) q' 
dK 3 + dK + dKK' + d 2 dq' 
qdd (V + dTK + ddK) 



Now, evaluating on the brane, we examine the order of 
each of these terms to find which is the dominant. For 
example, remembering that d and d y commute, we have 
for I > 0, 



4 l) ( dd2 l')\ y =o = d ° 2 (V n) (0)) ~ dK(0) 
and similarly 

B® (d (Yd + dT + T 2 ) q' + d K 3 + dK^ 

(dKK' + qdd (dTK + ddK) + 3 2 dq') 

Finally, 



y=o 
y=o 



dK(0) 
dK(0). 



3^(qT'3d) _ =dd J2 (?)g ( '- m) (0)r (i+1) (0), 



y=o 



and the I — term dominates this last sum, being of 
order K(0). Hence the dominant term in the expression 

for Ku is, as in the K„^ case, the one with the 
derivative of the Christoffel symbol, of the same order as 
K(0). 

We have now proved half of the inductive hypothesis, 
but still need to show that r( n+1 ) (0) - K(0). From (HJ|, 
we have 

r M Q „ (n+1) (y) = (d iC (n) ) + B„(d )-d a (d k$ ) 



-da ( , n) 



TK termsl 



(A17) 



rn 
n — 1 



dK^ m A 3^ (q a0 q w ) 



m=0 

By the inductive hypothesis 

dd^ (TK) ~ dK(Q), d^~ m) (qf p ) ~ dK(0) n>m 
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and 

d a (dK^ n) ) = d, a K£W + O(dK(0)). 

Therefore we can now read off from l|A17() the leading 
order contribution to T^ n \ 

r«J n+1 \0) = d,„K<ZW + d, v K«W-d' a K$ 

~ K(0), (A18) 

which agrees with l|A14jl and proves the inductive hy- 
pothesis. 

We are now finally in a position to calculate 
for general n. We know that the leading contribution 
obtained from repeatedly differentiating the right-hand 
side of l|A10p is 

KgW ~ q^T^-^dad 



which, from (|A18fl . immediately gives (|A2|) and the result 
is proved. 



for simplicity. The ij (with i ^ j) component of the Ein- 
stein equation, to first order in the perturbations, reduces 
to: 



$ - # - 2E' = ( 4— + 2% ) E, 

a a 



and the OUcomponent to: 



a 2d , 



(B3) 



(B4) 



So far these equations are equivalent to those one would 
have obtained in the low-energy limit. The difference 
comes from the 00-component of the perturbed Einstein 
equations: 

+ 1- (2f - df + d 3 g') $ - 2kH - 
a Ld V / 

k 2 d' 

—Ak 2 HE = a— cot d (2Sp + 3Sp) (B5) 

6 d 

and from the equation of motion for d: 



APPENDIX B: SCALAR PERTURBATIONS IN 
AN FRW BACKGROUND 

In this appendix, we shall present some of the details 
for the calculation of scalar perturbations of section lVBI 
with the metric perturbation as given in (|59|l . We recall 
that we picked the comoving gauge for which 5d — 0. In 
that gauge we then have: 

5\dd\ = d$ 

Terms that appear to be sub-dominant will only be 
dropped at the end. Using (j3*0|l . we get: 



SR = — (Sp - 3Sp) . 

lJ 



(Bl) 



Since a" = for the background (we assume the brane 
to be empty), (|Blll implies 



(B2) 



+ — (2k 2 £ + + 3*') 
a 

+ ^ (2* + 2E' - $) = ^-a 2 (Sp - 3Sp) 
3 oL 



We now perturb Ij31|l . writing 

z = \dd\ 
f(z) = ztanh(z/2) 



<!>' + (8^/ + | (/A/- -4/') )<!> (B(i) 



2k 2 E = — —a [cot d (28 p + 3Sp) - 3 coth d Sp 



Note that one must, at this order, treat and \E ,/ as 

four independent variables; differentiation with respect to 
conformal time will miss terms arising from higher order 
in ci, since d and d' are of different order. We must then 
solve the five equations (|B2IB0I simultaneously. Using 
(|B3(I to eliminate $ from l)B2j) . we obtain 

4fc 2 {^-n\E = 6 U 2fc 2 * + 6H ($' + 3*') (B7) 

k 2 

— —a 2 (Sp — 35p) 
±j 

We may use a combination of i|B5|) and (|B6|I to find an 
expression for $' in terms of <&,S!/,p and p and hence 
write E in terms of \&', p and p. This can then 
be used in i|B5|) to obtain a complicated expression for 
* in terms of 6 p and Sp. We then only keep the 

leading order in d for each coefficient, resulting in the 
much simplified equation 



d! 
d 



k 2 ^ = — a 2 - coth d Sp. 
6 d 



(B8) 



l(z) = - (tanh(z/2) +tan(z/2)) 

z 
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